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A shot noise, I(t), is a superposition of impulses occurring at random 
Poisson distributed times • • • , t-i , U , h , fe , • • • . In the simplest case, if 
the impulses all have the same shape F(t), then I(t) = ^2iF(t - U). We 
study, in this and more general cases, the distribution function Q(I) = 
Pr[I(t) ^ I]. One of our results is an integral equation for Q(I). This 
yields explicit expressions for Q(I) in a number of cases, including F(t) = 
c~ { ; it also permits a computational technique which is applied to F(t) = 
e~ sin ut for u> » 1. 

I. INTRODUCTION 

A shot noise, I(t), is a superposition of impulses occurring at random 
times • • • , t-i , to, k, t 2 , • • • . If the impulses all have the same shape, 
F(l), then 

I(t) = JlF(t - ti). (1) 

> 

More generally, the impulse shapes may be randomly chosen from a 
family of shapes, F(a, I), depending on a parameter a. Then 

lit) = £F(a,,* - U). (2) 

i 

We assume that the times t,- form a Poisson sequence with rate n impulses 
per second. In the case of (2), we assume that the parameters a, are 
chosen independently from a common distribution. 
We study the amplitude distribution function 

Q{I) = ErlJ(I) ^ /]. 

Rice 1 (Section 1.4) considered the noise ( 1 ) and noises (2) withF(a, t) = 
aF(t). He expressed the density function P(I) = Q'(I) as a Fourier 

333 
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integral. The Fourier integral is difficult to evaluate except by means of a 
series that Rice derived for the case of noises that are nearly gaussian 
(large impulse rate n). In our treatment we derive for (1) an integral 
equation 



di 



(3) 



f x dQ(x) - n { Q[I - F(t)]F(t) 

J— to J— oo 

or, equivalently, 

IQU) = f Q(x) dx + n f Q[I - F(t)]F(t) dt. (4) 

J — OO J— 00 

We solve (3) for some special choices of F(t) illustrated in Fig. 1. The 
section numbers on Fig. 1 give the part of the text in which each F (t) 
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Fig. 1 — Impulse functions F(l). 
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is considered. Analytic solutions are obtained for all cases except the last 
one of these; this case is an important one in practice, and so was chosen 
to illustrate a numerical solution of (3). 

For purposes of finding Q(I) any given noise of form (2) can be re- 
placed by an equivalent one of form (1). This equivalence is discussed 
in Section III. 

Several different ways of deriving (3) are possible. We give an analytic 
proof (Section IV) and a probabilistic proof (Section V). 

II. CHARACTERISTIC FUNCTIONS 

All impulse functions F(t), F(a, t) will be assumed integrable over 
— oo < / < oo . This assumption is no practical restriction and is made 
to ensure that the series (1) and (2) converge. 

We begin by deriving the characteristic function C(s) of 7(0 for the 
noise (2) ; i.e., 

C(s) = E[e-' H % (5) 

Here E denotes the expected value. 

The characteristic function C Q (s) in the case of the noise (1) is 

Co(s) = exp l-n j_ (1 - exp [-sF(t))) dt\ (6) 

(see Ref. 1, Section 1.4-7). We may obtain C(s) directly from (G) by 
regarding the noise (2) as a superposition of noises of the form (1) with 
different choices of F(i). Suppose, for example, that the parameter a 
has only a discrete range of possible values (a = Ai , A z , A 3 , • • • ), and 
let p k be the probability of picking a to be A k . Then, collecting to- 
gether the terms of (2) for which a, has the same value, one expresses 
I(t) as a sum of new independent random variables 

I(t) = h(t) +h(t) + ••• , 

where h(t) i« a noise of the form (1) in which the impulse shape is 
F(t) = F(A k , /) and the impulses arrive at an average rate, np k . per 
second. If (\-(s) is the characteristic function of h(t), 

C(s) = E[e-' wl)+ '* U)+ -\ 

C(s) = d(s)C,(s) ••• . (7) 

In (7) each C k -(s) may be evaluated by an expression of the form (6), 
and the final result is 

C(s) = exp [-/< (£ U " /^""""•"l} *)] • (8) 
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In (8) the expectation E is taken with respect to the random parameter 
a. Although our derivation of (8) used the assumption that a had a dis- 
crete range of values, a convincing limiting argument can be given for 
the truth of (8) in general. Alternatively one can rederive (8) in general 
by a slight modification of Rice's derivation of (6). 
In a similar way, we find 

E {exp [- sJ{ti) - •■■ - s n I(t n )]) 

= expl-n I 1 - E exp - jj s k F(a, t - r k ) dtV>, 
which might be used to study the joint distribution of /(ri), • • • , I(t n ). 

III. EQUIVALENCE 

In (6) it is evident that there are many different ways of choosing an 
F{t) to obtain the same distribution of I(t). The integral in (6) remains 
unchanged if F(t) is replaced by any other function F (t) such that, for 
every choice of U\ and u 2 , the two sets & and So of times t that are de- 
fined by: 

S: wi < Fit) ^ u 2 , 

(9) 
S Q : Mi < F Q (t) ^ u 2 

have the same measure: For example, the second function in Fig. 1 may 
be replaced by e -21 ' 1 . Some idea of the freedom with which one can con- 
struct such a new F (t) from a given F(t) may be had from Fig. 2. 
Given F(l) and n, one can construct a measure, dg(u), on the real u 
line by defining the measure of the interval ui < u ^ w 2 to be n times 
the Lebesgue measure of the set of times t for which u^ < F(t) ^ u. . 
Then, changing the variable of integration in (6) from t to u, one ob- 
tains 



CM = expl-f (1 - «—) <to(tt)]. 



(10) 





Fig. 2 — A pair of equivalent impulse functions. 
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Similarly, in the case of noise (2), let dg(u) be defined in such a way 
that the interval U\ < u ^ u* has measure equal to n times the expected 
value of the Lebesgue measure of the random set of times / satisfying 

ui < F(a, t) ^ u 2 . (11) 

Changing the variable of integration in (8) from t to u, one obtains for 
C(s) an expression that is just the right-hand side of (10). Thus, we 
call two noises which have the same dg(u) equivalent and have shown 

Theorem: The distributions of the amplitudes 7(0 of equivalent noises 
are the same. 

Given a noise (2) with measure dg(u), one can find a function F(t) 

such that 



a:*?*) 



/ ajrw) \ u forM>0 (12) 



and 



(-i:^) 



F - WW\ U ioru<0 . (13) 



The noise (1) with this choice of F(t) is equivalent to the given noise 
(2). Then every noise (2) is equivalent to a noise of form (1). For the 
problem of finding the amplitude distribution we now need consider 
only noises of the form (1). 

As a very simple application of the theorem, consider (2) with the 
family of impulse functions 

[6 if ^ t ^ a 

F(a,t) - | -. (14) 

[ otherwise, 

where the parameter a is distributed over positive values only. To find 
the measure dg(u), note that the set (11) has Lebesgue measure 



/.: 



dg(u) = 



oc if id < ^ Ui 

a if ^ Mi < 6 £t Ui 

otherwise. 



Then, dg(u) must lump all its measure onto two points u = and u = b. 
The measure of is oo and the measure of b is nA, where A = E(a), the 
expected length of the pulse (14). A noise of the form (1) that also has 
the measure dg( u) is the one with 

if ^ t £ A 

(15) 
otherwise. 
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This happens to be a noise for which the amplitude distribution is easily 
obtained. Since I(t) is just the number of impulses which arrive in the 
time interval from t — A to t, I{t) has the Poisson distribution 

Pr [I(t) = kb] = -^V"* e~" A . 

By the theorem, this result also solves the distribution problem for the 
original noise (14). In the same way, the example may be generalized 
as follows: If F(a, t) has the form S(t/a) for some given impulse function 
S(t), then the amplitude I(t) in (2) has the same distribution as the 
amplitude in (1) where F(t) is taken to be S[t/(E \ a | )]. 

Although the measure dg(a) determines the amplitude distribution, it 
does not determine all the statistical properties of the noise. One can 
easily find examples of functions F(t) and F (t) for which the corre- 
sponding noises (I), although equivalent, have different joint distribu- 
tions for the pair of random variables I(ti), 7(t 2 ). The spectrum of I(t) 
is proportional to the squared magnitude of the Fourier transform of F(t) 
(see Campbell's theorem in Ref. 1, Sections 1.2, 1.3), and so can be 
changed without changing dg(u). 

IV. DERIVATION 

A proof of (3) can be given from the formula (10) for the character- 
istic function C(s). We will assume here that the impulses F(l) or F(a, I) 
in question are nonnegative functions. This restriction will be removed 
in the next section and is made now in order to allow Laplace transform 
methods to be used. 

We have 

C(s) = E[e- 8H,) ] 



= f e-"dQ(I) 

= s f e-"Q{I) dl, 

Jo 



where the last formula is obtained by integrating by parts and noting 
Q(0—) = 0. Thus, C(s)/s is the Laplace transform of Q(I). 
Using (10), 

= — log C(s) 



C(s) ds 

= — I e~'"u dg(u). 
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Then 

_C\s) = CM( e -°» udg (u). (16) 

s s J 

Take the inverse Laplace transform of (16). The product on the right 
of (16) transforms into a convolution. One of the terms, C(s)/s, is the 
Laplace transform of Q(I). The other term is already in the form of a 
Laplace transform. Then the product transforms into 



/' 

Jo 



Q(I — u)u dg(u). 
This integral is another way of writing 

/ Q[I - F(t)]F(t) dt. 

To prove (3), it now remains to show that the term -C'(s)/s in 
(16) is the Laplace transform of 



•Ml 



x dQ(x). 



This follows because 

C'(s) = - f c~ SI I dQ(I) 

Jo 

and because 1/s is the Laplace transform of the unit step function. 

The integral equation (3) can be proved in several other ways; in 
particular, a more probabilistic proof will be obtained as a byproduct 
in the next section. 

V. SUMS OF NOISES 

Let Ii(t) and I t (t) be two independent shot noises with impulse re- 
sponses Fiit), F 2 (t), impulse rates n x , n 2 and measures dgi(u), dg 2 (u). 
Their sum I(t) = I,(t) + h{t) is a shot noise of the form (2). The im- 
pulse rate for the sum /(/) is n = n x + n 2 ; the function F(a, t) is 
chosen to be Fi(t) or F->(1) with probabilities ni/n and n ? /n. The meas- 
ure dg(u) for /(/) is 

dg(u) = dgi(u) + dg 2 (u). 
Since I(t) is the sum of the independent random variables h(t) and 
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h(t), the distribution of I(t) can be obtained from those of Ii(t) and 
^(0 by a convolution. 

This observation may be used occasionally to compute the distribution 
function for a given noise I(t). Suppose the measure dg(u) of I(t) can 
be written as a sum 2 dgi(u) of measures of noises /,•(£) for which the 
distribution functions are known. Then 1(1) is equivalent to a sum of 
independent noises 2 Ii(t) and its distribution function may be ob- 
tained by convoluting the distributions of h(l) together. Even when an 
exact decomposition dg(u) = 2 dgt(u) is not known one might approxi- 
mate dg(ii) by a sum to get an approximate Q(I). 

A particular instance of a decomposition is the following. Let meas- 
ures dg (u) and dg~(u) be defined by 



dg + (u) = 
dg~(u) = 



dg(u) if u > 

otherwise, 

^dg(u) if u ^ 

otherwise. 



Then dg(u) = dg + (u) + dg~(u). At present, our proof of (3) holds 
only for positive noises. A similar derivation of (3) for negative noises 
also holds. In the general case we might consider our noise I(t) to be 
equivalent to a sum I + (l) + I~{t) of independent positive and negative 
noises, and compute Q(I) as a convolution 



QU) = / Q + U - x) dQ-(x). 



That (3) holds in general now follows from the next lemma. 
Lemma: Let dgi(u), dg*(u), dg{u) be measures such that 

dg(u) = dgi(u) + dg 2 (u). 

Let Qi(I), Qi(I) be solutions of the integral equation corresponding to the 
measures dgi(u), dg*(u). Then for measure dg(u) the convolution 

QU) = / Qi(/ - w) dQ 2 (w) 

= / Qt(I - x) dQM 
is a solution of the integral equation (3). 
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Proof: 

( y dQ(y) = [ f y dQ x {y - w) dQ 2 (w) 

J— oo J—x J w 

= [ f [(y - w) + w] dQi(y - w) dQ 2 (w) 

J— oo ■'ui 

= / / Qi(I — w — u)u dgi(u) dQ 2 (w) 

J— oo w XV 

+ / / dQi(y - w)Q 2 (I - w)u dg 2 (u) 
= f Q(/- m)m^(m). 

J— 00 

Taking rffifi and rf^» as rfgf + and cty - , the lemma, together with the result 
of Section IV, completes the derivation of (3). 

A different proof of (3) may now be outlined as follows. Consider 
first the noise with response function (15). The integral equation (3) 
holds in this special case. For dg(u) gives measure nA to the point u = b 
and 

QU) = £ ^.-'. 

kb § / ft- I 

Then 

J-oo kb£l Kl 

= nAbQ(I - b) 

= / Q(I - u)u dg(u). 

J— oo 

Next, any step function F(t) with a finite number of steps has a meas- 
ure dg(u) = 2 dgi(u), where each dg%(u) concentrates its measure on a 
single value of u [a level of F(t)]. We have just proved that the inte- 
gral equation holds for each of the dgi(u) noises. By the lemma we can 
conclude that (3) holds for all step function noises. By limiting argu- 
ments, one might establish (3) more generally. 

VI. EXAMPLES 

For certain choices of F(l) the integral equation can be solved easily. 
Some special cases of this kind will be examined in this section. 
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6.1 Example 1 

First consider a noise ( 1 ) with 

(e~ l if t £ 
F(t) = 

[0 if t < 0. 

We expect a distribution Q(I) that has a density P(I) = Q'(I) and 
so write the integral equation in the form 

IP(I) = n [ P (I - u) du = n [ P(x) dx. 

Differentiating, we obtain 

IP'(I) - (n - l)P(I) = -nP(I - 1). 

To solve this differential difference equation, note first that, when 
0|/<1,P(/-1)=0. Hence, for ^ / < 1, 

where c is a constant of integration to be determined. For larger values 
of /, the differential difference equation may be converted to an integral 
form 

p(I) = 7'- 1 \c-nj P(x - l)x~" dx\ . 

Since the integrand is known for x < 2, we can determine P(I) for 
I < 2. Next, this result enables us to integrate further to get P(I) for 
I < 3, etc. Clearly, the analytic form of P(I) changes at / = 1, 2, 3, 
• ■ •. For example, when n = 1, we have 

P(I) = c ifO £ / £ 1, 

P{I) = c(l - log/) if 1 ^ i" g 2, 

P{I) =c ri-log/ + £ log( ^- 1) ^] if 2 ^ J ^ 3, 

etc. 
Finally, the constant c must be determined by the condition 

[ P(I) dl = 1. 
Jo 

The constant can be determined as follows: The Laplace transform 
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pis) of P(I) is the characteristic function. By (6), 

pis) =exp(-n/ Q l -^dy). 
This may be rewritten, with the aid of partial integration, as 
■p{s) = exp — n(l — e") log s + 

n \ e~ y log y dy — n I e~" log y dy 

Jo Jj 

= S - n e- ny {l + 0[r 8(1 "°]} for any e > 0. 
Thus, for < I < I, 

w > - m r "' 

where 7 = 0.577215665- • • is Euler's constant. Hence 

c = 



T(n)' 

6.2 Example 2 

Our next example concerns a noise (1) with Fit) defined as follows: 
Fit) = for t ^ 0, 
and 

,.r e* !or i>o. 

•>f(» y 

This somewhat artificial noise interests us because it has a very simple 
F(7). A sketch of Fit) is shown in Fig. 1, marked "Section 6.2". For 
small t, Fit) grows large but only logarithmically: 

F{t) ~ -log t. 

For large /, F(t) has an exponential tail, 

Fit) ~y\ 

where 7 is Euler's constant. 

For this noise, dgiu) = ne~ u du and the integral equation can be put 
in the form 



I Pine 1 = n [ Pix)e x dx. 

Jo 
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Differentiating, we obtain a very simple differential equation and the 
solution is 

-Ijn-l 
PU) = € 



T(n) 

This solution exhibits a rapid approach of P(I) to gaussian form (with 
mean n and variance n) as n — * °° . 

0.3 Example 3 

If we consider the case 

f 1 - t ^ t ^ 1 

F(t) = \ 

[ otherwise, 

then the distribution Q(I) has a jump of e~ n at / = 0, since e~ n is the 
probability that no point of the Poisson process falls into an interval of 
unit length. We therefore seek a density function P(I) such that 

Q(I) = e ~ n + [ P(x) dx for 7^0. 
Jo 

The integral equation then becomes 

n \ P(I - y)y dy + nle~" if < I g 1 
Jo 

n [ P(I - y)y dy if 1 < I. 



mi) = 



'0 

If P(I) = R"{I), this becomes 

[0 if/<l 

IR (I) - nR(I) = 

I -n[R'(I - 1) + R(I - 1)] if I £ 1. 

This can be solved recursively. In the first interval (I < I) 

where 7i is the Bessel function, and where the coefficient of I x has been 
determined by substitution in the integral equation. 
In the general case, if I > 1, 

P(/) = e -«£^V +1),2 (/ - fc) a+1,,2 / i - 1 [2Vn(7 - *)], 
k=a K ! 
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where [I] is the largest integer ^/. For example, if 1 < / ^ 2, 
P(I) = <T" [/j/f /i(2Vn/) " n/o(2Vn(/ - 1))J . 

(Note that I-\ = /i .) These formulas can be derived either from the 
integral equation, or, with rather more courage but actually less work, 
from (6). 

6.4 Example 4 

Choosing 

f -logt, < I ^ 1 
F(l) = 

[ otherwise 

provides another simple case. Again, Pr(7 = 0) = e~", and we seek a 
density P(I), as in Example 3. For values I > 0, P(I) satisfies 

IP(I) = n f P(I - u)ue~" du + ne~ n Ie~ l . 

.'o 

Now, letting R(I) = Ie'P(I), we obtain again the differential equation 
IR"{I) + nR(I) = 0, this time for all I > 0. The .solution is 

P(I) = <f </+ " ) /j/f /i(2\/^) for / > 0. 

6.5 Example 5 

In this example we let F(t) be the response of a simple tuned circuit 
to an impulse; i.e. 

F(t) =e-'sm(*f). (17) 

Although the period h appears in F(t) as a parameter, the correspond- 
ing measure dg h (u) tends to a limiting measure dg{u) as h — > 0. We will 
solve the integral equation numerically in this limiting case only. We 
then expect this result to be applicable as a good approximation when- 
ever the tuned circuit has high Q; i.e., h « 1. 

To get the limiting measure dg(u) let us examine F(t) in a small 
neighborhood of a time t = T at which the sinusoid is at a maximum. 
During the period from T to T + h the exponential e~' changes only 
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by a factor e~ h = 1 + 0(h). Thus, aside from terms of order 0(h), 
F(t) is 

F(t) = e cos 

in this period. Given a level u > 0, F(t) ^ u for a time 

- arc cos (ue r )[l + 0(h)] 



IT 



during the period T 2a t ^ T + h. For small h, we conclude that F(l) 
lies above level u for a total amount of time approaching 

Jl) W 7T Jo 

1 /* arc cos Z irr 

= ~ / t rf^ 

IT J u % 

as /i — > 0. Similarly, setting the amplitude level u at a negative value, 
we obtain 



f u dg(w) _ 1 f arc cos Z ,„ 

J-oc W 7T J|«| 2/ 



The measure dg'Cw) is now known. We wish to solve the integral equa- 
tion in which 



dg(u) = 



f n arc cos I u I rf?t ... , „ . 

■" I M I (18) 

if | u | ^ 1. 



In this case, we had to resort to numerical methods. The integral 
equation might be approximated directly by a system of linear algebraic 
equations. However, such an approximation would be troublesome in 
our case because the integral equation is homogeneous. Unless we could 
guarantee that the approximating system would have determinant zero, 
there would be no nontrivial solution at all. The procedure that fol- 
lows avoids this difficulty. 

Let measures dg + (u), dg~(u) be defined as in Section V. We will solve 
the two integral equations with measures dg (u) and dg~(u) and con- 
volute the two solutions P + (I), P~(I) together to get the desired den- 
sity P(I). This approach has the advantage that the integral equation 
expresses P + (I) in terms of only values of P + for arguments <I. Thus, 
one can compute P + (I) approximately by a simple recurrence. 



AMPLITUDE DISTRIBUTION OF SHOT NOISE 



347 



One might start the recurrence computation by assuming a value of 
P + (I) for a small /; afterward the solution could be normalized to make 



f P + (I) dl = 1. 

Jo 



However, we can obtain in the Appendix an asymptotic formula for 
P + (I) for / near 0, and thereby start the recurrence with nearly correct 
values of P + (I). 

Figs. 3 and 4 give the results of the computation for this limiting case 
of an infinitely rapidly oscillating tuned-circuit response. Computations 
were made for rates of 2.5, 5 and 10 impulses per unit time, where the 
time scale is determined by the exponential in (17). Fig. 3 shows P(I) 
for these three cases; Fig. 4 plots Q(I) on log normal paper and com- 
pares the result with the gaussian of the same mean and variance. As 
expected from Rice's theory, the noticeable differences are in the tails 
of the distributions; in the shot noise, very large values of / are more 
likely than in the corresponding limiting gaussian. 
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Fig. 3 — Amplitude density P(I) for high-frequency damped sinusoid noise; 
F(l) = e~< sin at, u > 1; P(-'l) = P(D- 
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Fig. 4 — Comparison of amplitude distribution Q{I) with gaussian approxi- 
mation; Q(-I) = 1 - <?(/). 

As indicated above, the computation was performed in two stages. 
First, the integral equation for P + (I), 



IP + (I) = f P + (I - u) dg(u), 
Jo 



was solved by approximating the integral through Simpson's formula 
and solving the resulting triangular system of equations; the series for 
P + (I) near 1 = was used to start the computation. Then, with the 
aid of the theory in Section V, P(I) was computed as 

P(I) = P P + (y)P + (I + y) dy, 

and Q(I) as the integral of P. 



VII. ACKNOWLEDGMENT 



We wish, at this point, to thank Miss M. A. Lounsberry, for writing 
a general program for the IBM 704 for solving integral equations of 
this kind. Without her ingenuity and perseverance the present numerical 
results would not have been obtained. 



AMPLITUDE DISTRIBUTION OF SHOT NOISE 349 



APPENDIX 



Derivation of the Form for P (I) for Small I for a Tuned-Circuit Response 
We begin with the integral equation 

i>min(l,7) 

IP + (I) =- P + (I - u) arc cos udu. 

IT Jo 

Taking Laplace transforms, we obtain 

p(s) +|^l[l - Us) 4- Us)} = 0, 

where 7„ is the modified Bessel function, and L the Struve function. 
Remembering the condition p(0) = 1, and that 7 (0) = 1, L (0) = 0, 
we find 

>U\ T n f 1 - Zo(tt) + Lq(m) ^„"| 

Now 7o — 7/o approaches zero at infinity rapidly enough for that portion 
of the integral to converge as s — > » , but the remaining term diverges. 
We thus write 

r .. /•• r /..\ r /.^ i („. i i \-n 



w x r n r l «( m ) - 7 °(«) + (« +_i) n 

p(«) = exp|_-- J o J 

«p[» JT U.Q-^O + U + i)- 1 dtt j 



•exp^, 
By Ref. 2, p. 426, 



7 (u) - L (m) 2 f x J (x) 



Wo a- 2 + w 2 



u 

Hence, 

7 (m) - Lo(tt) - (m + l)" 1 2 f/o(*) - (1 + X 2 )" 1 



_ 2 r 

7T ^0 



dx. 



If T Jo X 2 + tt 2 

We may now integrate with respect to u and interchange integrations. 



350 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH I960 



There results 

'•/ooo-iooo-u+ir 1 . rjo(x) - a + x 2 y 



r h(u)-uu)-(u + ir l du = r 

Jo it Jo 



dx 



I [ Ux) ~ (TtWJ log ■ * 



= log 2-7 
by Grobner and Hofreiter' and Bierens de Haan. Thus, 
p(s) = C - ( " /2)T (2)" /2 ( S + l)-" /2 

\n r U(u) -7p(u) + (i<+ I)" 1 1 

It remains only to discuss the behavior of the last integral for large s. 
By Ref. 2, p. 332, 



2 /"* sin ux 
Lo(w) - h(u) = — - I „/ , , , dx. 

T Ja V 1 + X 2 



Hence, 



while 



J« u t Jo VI + x % J sx y 

= - / log(z + VI + s 2 ) ■ dx, 

IT Jo X 

I. ^rh) = In ( 1+i s )' 



Both of these are 0(l/s) for large s, and hence, for large s, 

p(s) = e-'^-s-" 12 + 0[s- M2) -\ 
Hence 

p+(I) = 2e (21) + 0(/n/2) 



(1) 



for small values of I. 
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